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Electron transport experiments on two lateral quantum dots coupled in series are reviewed. An 
introduction to the charge stability diagram is given in terms of the electrochemical potentials 
of both dots. Resonant tunneling experiments show that the double dot geometry allows for an 
accurate determination of the intrinsic lifetime of discrete energy states in quantum dots. The 
evolution of discrete energy levels in magnetic field is studied. The resolution allows to resolve 
avoided crossings in the spectrum of a quantum dot. With microwave spectroscopy it is possible to 
probe the transition from ionic bonding (for weak inter-dot tunnel coupling) to covalent bonding 
(for strong inter-dot tunnel coupling) in a double dot artificial molecule. This review on the 
present experimental status of double quantum dot studies is motivated by their relevance for 
realizing solid state quantum bits. 
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I. INTRODUCTION 



Quantum dots are man-ma de sub-micron structures i n a solid, typically consisting of 10 3 -10 9 atoms and a com- 
parable number of electrons (Kouwenhovcn et al, 1997). In semiconductor quantum dots all electrons are tightly 
bound, except for a small number of free electrons, which can range from zero to several thousands. For the quantum 
dot devices considered in this review, the starting point for fabrication is formed by a heterostructure consisting of 
different semiconductor materials (GaAs/AlGaAs). The free electrons are strongly confined to the interface between 
GaAs and AlGaAs, forming a 2-dimensional electron gas (2DEG). Confinement in the other two dimensions is ac- 
complished by locally depleting the 2DEG, via etching techniques or metal gate electrodes. The resulting structure 
is weakly coupled to source and drain electrical contacts by tunnel barriers, connecting it to the outside world. 



An important clement of electronic transport through quantum dots is Coulomb blockade ( Averin and Likharev^ 
198<j|, 1991; Grabert and Devoret, 1992; Kouwenhovcn et al, 1997). An extra electron can only be added to the dot, 



if enough energy is provided to overcome the Coulomb repulsion between the electrons. Next to this purely classical 
effect, the confinement in all three directions leads to quantum effects that strongly influence electronic transport at 
low temperature. In particular it leads to the formation of a discrete (OD) energy spectrum, resem bling that of a n 
atom. This and other similarities have therefore lead to the name artificial atoms for quantum dots (Kastner, 1993). 

The next logical step after studying individual quantum dots is to study systems of more than one dot. Where single 
quantum dots are regarded as 'artificial atoms', two quantum dots can be coupled to form an 'artificial molecule'. 
Depending on the strength of the inter-dot coupling, the two dots can form ionic-like (weak tunnel coupling) or 
covalent-like bonds (strong tunnel coupling) . In the case of ionic bonding the electrons are localized on the individual 
dots. The binding occurs, because a static redistribution of electrons leads to an attractive Coulomb interaction. 
Weakly, electrostatically co upled quantum dots with neg ligible inter-dot tunnel conductance are covered by orthodox 
Coulomb blockade theory (Averin and Likharev, 1991). In the case of covalent bonding, two electron states are 
quantum-mcchanically coupled. The main requirement for covalent binding is that an electron can tunnel many times 
between the two dots in a phase-coherent way. Here the electron cannot be regarded as a particle that resides in one 
particular dot, but it must be thought of as a coherent wave that is delocalized over the two dots. The bonding state 
of a strongly coupled artificial molecule has a lower energy than the energies of the original states of the individual 
dots. This energy gain forms the binding force between the two dots. 

The theoretical possibility to perform certain tasks in a much more efficient way using a 'quantum computer' instead 
of a 'classical computer', has stimulated the search for physical realizations of the basic building block of such a com- 
puter: the quantum bit. In principle, any quantum two-level system can be used as such a qubit. In particular, recent 



studies have put forward double quantum dots as interesting candidates for realizing qubits (Loss and DiVincenzo 
199c|). The possible application of double quantum dot devices in quantum logic forms an important motivation for 
this work. 

In this review we concentrate on electron transport through lateral double quantum dots coupled in series. All 
devices have been fabricated and all experiments have been performed at Delft University of Technology and NTT 
Basic Research Laboratories. By now there exists an extensive literature on experimental studies of electron transport 
through double lateral quantum dots coupled in series 1 , and lateral double dots coupled in parallel 2 . Vertical double 
quantum dot structures 3 fall outside the focus of this review. In vertical structures, the characteristics of the tunnel 
barriers are set by the growth parameters of the heterostructure, limiting the experimental tunability. Besides that, 
the gate geometry used in these devices, makes it difficult to address dots independently. 

As a first step to understanding double dot systems we introduce the stability diagram (Pothicr et al, 1992), or 
honeycomb diagram, in section 0. It is a convenient tool in the analysis of double dot transport properties. Resonant 
tunneling experiments discussed in section [II, show that the resonant widths are only determined by the lifetime of 
the discrete energy states, independent of the electron temperature. In section IV we discuss level spectroscopy in 
a magnetic field. The double dot geometry offers sufficient energy resolution to probe intra-dot level repulsion. In 
section [v| we present microwave spectroscopy measurements on a quantum dot molecule. We illustrate the transi- 
tion from a weakly coupled double dot to a strongly coupled double dot, by discussing a two-level system in section 
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V.A. Although being a clear simplification, the mapping of the double dot on a two-level system grasps much of the 
physics of the experiments presented in section Irradiation with microwaves leads to photon assisted tunneling 
(PAT) (sections V.B and V.C), which turns out to be a powerful tool not only to reveal the character of the inter-dot 
coupling, but also to quantitatively determine the bonding strength. 



II. STABILITY DIAGRAM 



In this section we introduce the stability - or honeycomb - diagram that visualizes the equilibrium charge states of 
two serially coupled dots. 



A. Linear transport regime 



1. Classical theory 



We start with a purely classical description i n which the influen ce of discrete quantum states is not taken into 
account yet (Dixon, 1998; Pothier et aL, 1992 ; Ruzin et al, 199~2] ) 4 . The double dot is modeled as a network of 
resistors and capacitors (Fig. |l|). The number of electrons on dot 1(2) is Ni^)- Each dot is capacitively coupled to 




FIG. 1 Network of resistors and capacitors representing two quantum dots coupled in series. The different elements are 
explained in the text. Note that tunnel barriers are characterized by a resistor and a capacitor, as indicated in the inset. 



a gate voltage V g i^) through a capacitor C s i(2) and to the source (S) or drain (D) contact through a tunnel barrier 
represented by a resistor Rur) and a capacitor Cl<r\ connected in parallel. The dots are coupled to each other by 
a tunnel barrier represented by a resistor R m and a capacitor C m in parallel. The bias voltage, V, is applied to the 
source contact with the drain contact grounded (asymmetric bias). In this section we consider the linear transport 
regime, i.e. V ~ 0. If cross-capacitances (such as between V g \ and dot 2), other voltage sources and stray capacitances 
are negligible, the double dot electrostatic energy reads (a full derivation is given in the appendix) 

U{N X ,N 2 ) = \NlE cl + X -NlE C2 + N^Ecm + f(V gl , V g2 ) (1) 
f(V gl ,V g2 ) = ^-{C gl V gl {N 1 E cl + N 2 E Cm ) + Cg 2 V g2 (N 1 E Cm + N 2 E C2 )} 

+ ~^{~^ c2 g\Vg\ E c\ + -jC 2 g2 V g2 Ec2 + C gl V gl Cg 2 Vg 2 Ec m } 



where £'ci(2) is the charging energy of the individual dot 1(2) and Ec m is the electrostatic coupling energy. The 
coupling energy Ec m is the change in the energy of one dot when an electron is added to the other dot. These energies 
can be expressed in terms of the capacitances as follows 



4 Ruzin et al. have an eror in their Eq. (9), namely the numerator, (2C + CiSij). For i = j, this relation gives (2C + Ci) when it should 
read (2C + C2), and vice- versa. 



Here Cu2) is the sum of all capacitances attached to dot 1(2) including C m : Cu 2 ) — ^L(B) + Cgi(2) + C m . Note that 
Eqi(2) can be interpreted as the charging energy of the single, uncoupled dot 1(2) multiplied by a correction factor 
that accounts for the coupling. When C m = 0, and hence Ec m — 0, Eq. (Q) reduces to 

V(N U N 2) = + <W>' + l-*M+C*V*r (3) 

zCi 2C 2 



This is the sum of the energies of two independent dots. In the case when C m becomes the dominant capacitance 
(C m /C 1 ( 2 ) — » 1), the electrostatic energy is given by 

U(N U N 2 ) = [-(^ + ^\+C^ + CM 2 (4) 

2(Ci + C 2 ) 



This is the energy of a single dot with a charge iVi + N 2 and a capacitance of C\ + C 2 , where Cu 2 ) = Ci(2) — Cm is 
the capacitance of dot 1(2) to the outside world. Thus, a large inter-dot capacitance C m effectively leads to one big 
dot. 

The electrochemical potential £ti( 2 )(iVi, N 2 ) of dot 1(2) is defined as the energy needed to add the iV^jth electron to 
dot 1(2), while having -/V 2 (i) electrons on dot 2(1). Using the expression for the total energy Eq. (|l]), the electrochemical 
potentials of the two dots are 



^(Nx, N 2 ) = U(N U N 2 ) ~ U(N X ~ 1, N 2 ) 

= (JVi - ~)£ C i + N 2 E Cm - {C gl V gl E C i + C g2 V g2 E C m) 



(5) 



ia{Nx,N 2 ) = U(N U N 2 ) - U{N U N 2 - 1) 

= (N 2 - X -)E C2 + NxEcm ~ t~t {CgiVgiEcm + C g2 V g2 E C2 ) 



(G) 



The change in fi 1 (Ni,N 2 ) if, at fixed gate voltages, N-y is changed by one, /ii(iVi + 1, N 2 ) — [j,i(Ni,N 2 ) — Eci, is 
called the addition energy of dot 1 and equals the charging energy of dot 1 in this classical regime. Similarly, the 
addition energy of dot 2 equals E C2 , and fx^Nt, N 2 + 1) - (J,i(Ni,N 2 ) = fJ, 2 (Nt + 1, N 2 ) - fi 2 (N l7 N 2 ) = E Cm . In 
the next section we will discuss the addition energy in the quantum regime, where also the spacing between discrete 
energy levels plays a role. 

From the electrochemical potentials in Eqs. (||) and (||) we construct a charge stability diagram, giving the equilib- 
rium electron numbers iVi and N 2 as a function of V g i and V g2 . We define the electrochemical potentials of the left 
and right leads to be zero if no bias voltage is applied, /i^ = fiji = 0. Hence, the equilibrium charges on the dots are 
the largest values of A^i and N 2 for which both /zi(iVi, N 2 ) 

and fj, 2 (Ni, N 2 ) are less than zero. If either is larger than zero, electrons escape to the leads. This constraint, plus the 
fact that iVi and N 2 must be integers, creates hexagonal domains in the (V g i, V^-phase space in which the charge 
configuration is stable. 

For completely decoupled dots (C m = 0) the diagram looks as in Fig. ||(a). The gate voltage V^i( 2 ) changes the 
charge on dot 1(2), without affecting the charge on the other. If the coupling is increased, the domains become 
hexagonal (Fig. gj(lo) ) . The vertices of the square domains have separated into 'triple-points'. When C m becomes the 
dominant capacitance (C m /Ci( 2 ) — > 1), the triple-point separation reaches its maximum (see Fig. ||(c)). The double 
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FIG. 2 Schematic stability diagram of the double-dot system for (a) small, (b) intermediate, and (c) large inter-dot coupling. 
The equilibrium charge on each dot in each domain is denoted by (Ni, N2). The two kinds of triple-points corresponding with 
the electron transfer process (•) and the hole transfer process (O) are illustrated in (d). The region in the dotted square in 
(b) is depicted in more detail in Fig. [5]. 



dot behaves like one dot with charge Ni + N 2 , as seen from Eq. (||). 

We are considering the linear regime of conductance, implying /i^ — /i/? = — |e|V « 0. In order to obtain a measurable 
current, the tunnel barriers need to be sufficiently transparent. At the same time, however, the tunnel barriers need 
to be sufficiently opaque to ensure a well-defined electron number on each dot. For double dots coupled in series, a 
conductance resonance is found when electrons can tunnel through both dots. This condition is met whenever three 
charge states become degenerate, i.e. whenever three boundaries in the honeycomb diagram meet in one point. In Fig. 
||(d) two kinds of such triple-points are distinguished, (•) and (O), corresponding to different charge transfer processes. 
At the triple point (•), the dots cycle through the sequence (Ni,N 2 ) -» (Ni + 1,N 2 ) -» (Ni,N 2 + 1) -» (Ni,N 2 ), 
which shuttles one electron through the system. This process is illustrated by the counterclockwise path e and the 
diagram of an electron sequentially tunneling from the left lead to the right in Fig. ||(d). At the other triple-point (O), 
the sequence is (Ni + 1, N 2 + 1) -> (Ni + 1, N 2 ) -> (A^i, N 2 + 1) -> (iVi + 1, N 2 + 1), corresponding to the clockwise 
path h in Fig. @(d). This can be interpreted as the sequential tunneling of a hole in the direction opposite to the 
electron. The energy difference between both processes determines the separation between the triple-points (•) and 
(O), and is given by Ecm, as defined in Eq. (|). 

The dimensions of the honeycomb cell (see Fig. ^) can be related to the capacitances using Eqs. (g) and (g). From 

fx 1 (N 1 ,N 2 ;V gl ,V g2 )=fx 1 (N 1 + l,N 2 ;V gl + AV gl ,V g2 ) (7) 

we obtain 

*V gl = P- (8) 

and similarly we can derive 

W g2 = P- f ) 
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AV g2 m 









FIG. 3 Schematic stability diagram showing the Coulomb peak spacings given in Eqs. (||) and ([HJ). These spacings can be 
determined experimentally by connecting the triple-points. 



From 



we obtain 



»i(N 1 ,Nr,V gl ,V g2 )=(Mi{N 1 ,N 2 + l;V gl +AV g v ?, V g2 ) 



AT/TO _ M^™ 



= AKi 



C m 



and similarly we can derive 



AV, 



\C„ 



<)2 



C g2 C\ 



AV, 



y2 



Ci 



(9) 



(10) 



However, for a full characterizatio n of a ll capacitances in the system an analysis in the non-linear transport regime is 



required, as is discussed in section II. B 



Above we assumed that V g \ and V g2 only couple directly to the respective dots. In practice, however, there is a 
finite cross-capacitance from one gate to the other. The respective cross-capacitances result in a change of the slope 
of the charge domain boundaries in the honeycomb diagram. From Figs. g(b) and || it is clear that both kinds of 
triple points (• and O) form a square lattice. However, with finite cross-capacitances the positions of the triple points 
move to lower V g i($\ for increasing V g2 m (at constant V). 



2. Quantized states 

The discussion of the stability diagram so far has been completely classical. However, the strong confinement of 
electrons in the dots can lead to the formation of a discrete energy spectrum. To account for the quantized energy 
states in the dot, we need to incorporate their energies in the electrochemical potential. The electrochemical potential 
for adding an electron into energy level n of dot i is denoted by Hi <n . Within the constant interaction model, Hi <n is 
the sum of the classical electrochemical potential fj,f ass and the single-particle energy E n : fj,i jn = fj l ( r lass -\- E n . In the 
classical regime we found that the addition energy (the change in electrochemical potential needed to add an extra 
electron) equals the charging energy Eci (for dot 1) or Eq 2 (for dot 2). In the quantum regime, the addition en- 
ergy for the (A^i + l)th electron occupying discrete level m, with the iVith electron occupying discrete level n, becomes 



(Ni + l,N 2 ) - m,n(Ni,N 2 ) = Eci + (E m — E n ) 

= Eci + AE. 



(11) 
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Similarly, we find Eq2 + AE for the addition energy of dot 2. Note that for a (spin-) degenerate level AE can be zero. 
The dimensions of the honeycomb cell as given in Eqs. (||) and (|l0|), and depicted in Fig. |3] for the classical regime, 
change as follows 



< 



AE 



AV gm =-^(l + - (12) 



C ffl (2)C 2 (i) V ECm 



The electronic configuration that gives the lowest possible total energy in dot 1(2), is referred to as the dot 1(2) 
ground state. Any configuration with a higher total energy is referred to as an excited state. The electrochemical 
potential for adding the A^ 1 ( 2 )th electron to the lowest unfilled energy level of the (Nu2) — l)-electron ground state is 
labeled fiifl(N±, N2) {a*2,o(-^1) -^2)}. The electrochemical potential for adding the A^ 1( - 2 )th electron to a higher unfilled 
level of the (iVi(2) — l)-electron ground state - or to any unfilled level of an (N1/2) — l)-electron excited state - is 
labeled ^(JV^JVa), Mi, 2(^1,-^2), - {^2,1(^1,-^2), ^2,2(^1,^2), ••■}• 




FIG. 4 Schematic diagram of the electrochemical potentials fn, n {Nx, N%) in dots and leads in the linear regime. The first 
subscript indicates either the lead (L,R) or the dot (1,2). The second subscript refers to the nature of the dot energy state 
(ground state, n — or nth excited state. 



In Fig. |J a schematic diagram is given, showing the electrochemical potentials in the leads and dots in the linear 
regime — Ma. — ~ 0)- The ground state electrochemical potentials ^1^(1,0) and /i2,o(0, 1) align within 

the small bias window, allowing an electron to tunnel from left to right. This is an example of an electron transfer 
process as depicted in Fig. |^(d) . Note that an alignment of an arbitrary combination of electrochemical potentials in 
dot 1 and dot 2 does not necessarily lead to a current. For example, the alignment of Mi,o(l, 0) and M2,o(l, 1) does not 
result in current through the double dot. In the linear regime electron transport occurs via ground states, whereas 



the excited states start to play a role in non-linear transport, as will be discussed in section [LB. In the following 



discussion of the linear regime the ground state of dot 1(2) is denoted by Mi(2) (Ni, N2) (without the discrete level 
index) . 

A more detailed picture of a honeycomb cell in the linear regime, marked by the dashed square in Fig. |2|(b), is given 
in Fig. ||. The configuration of the ground-state electrochemical potentials is given in schematic diagrams on some 
places in the stability diagram. The dashed lines, which are extensions of the solid lines forming the honeycomb cells, 
help to find the position of both electrochemical potentials on a certain place. A crossing of dashed lines (as in the 
charge domains (0,1) and (1,0)) indicates that two electrochemical potentials align, but does not result in a current 
through the double dot. 
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FIG. 5 Region within the dotted square of Fig. |2|(b), corresponding with the 'unit cell' of the double-dot stability diagram. 
Four different charge states can be distinguished, separated by solid lines. At the solid line connecting the two triple points, 
the charge states (0,1) and (1,0) are degenerate. At the other solid lines the electrochemical potential of at least one dot is 
zero and thus equals the electrochemical potential of the leads. The dashed lines are the extensions of the solid lines within 
the honeycomb cells. The triple-points lie on the crossing points between the solid lines. The schematic diagrams show the 
configuration of the ground-state electrochemical potentials on the corresponding place in the honeycomb diagram. 



3. Experimental stability diagrams 

Before discussing some experimental stability diagrams, we introduce the two kinds of lateral double dot devices 
being studied in this review. For the first type, only metal gate electrodes are used to confine the electrons in the 
2DEG beneath. For the second type we use a combination of metal gates and dry etching to realize confinement. 

A scanning electron microscope (SEM) image of the first device is shown in F ig. |](a). Metal gates are d eposited on 



top of a GaAs/AlGaAs heterostructure with a 2DEG 100 nm below the surface ( [van der Vaart et ai, 1995|) . Applying 
a negative voltage to all gates depletes the 2DEG underneath them and forms two quantum dots. Current can flow 
from the large electron reservoir on the left via the three tunnel barriers induced by the gate pairs 1-F, 2-F, and 3-F 
to the reservoir on the right. The transmission of each tunnel barrier can be controlled individually by the voltage on 
gates 1, 2 or 3. A single quantum dot can be defined in the 2DEG by applying only a voltage to gates 1, 2, I and F 
(dot 1) or to gates 2, 3, II and F (dot 2). In this way, the individual dots can be characterized and their properties 
compared to those of the double dot. 

The second device is schematically shown in Fig. ^](b). First a channel is defined in the 2DEG by focused-ion-beam 
( FIB) or electron-cyclotron-r esonance (ECR) etching of an Alo.3Gao.7As/GaAs modulation-doped heterostructure 



(Fujisawa and Tarucha, 1996). A double quantum dot can be formed by applying negative voltages to gates GL, GC 
and GR. 

All experiments have been performed in a dilution refrigerator with a base temperature of 10 mK. The effective 
electron temperature in the leads is higher and can vary between ~40 and ~100 mK. A significant source of heating 
is the noise coming from the measurement electronics. The filters, used to attenuate the noise, have to be effective 
over a very large band width. They consist of a distributed RC network, usually a thin resistive wire going through a 
conducting medium such as copper powder or silver epoxy. The filters are installed at low temperature to minimize 
the thermal noise of the resistors inside. The filters are integrated with the sample holder in such a way that all 
sample wires are carefully shielded once they are filtered. 

To effectively create a double quantum dot, all gate voltages need to be tuned properly. The cross capacitances 
between the various gate electrodes, make it difficult to vary just a single parameter without affecting the others. 
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FIG. 6 Double quantum dot devices, (a) SEM micrograph of a double dot defined by metallic gates (light gray areas). The 
ungated 2DEG (100 nm below the surface) has a mobility of 2.3 x 10 6 cm 2 /(Vs) and an electron density of 1.9 x 10 15 m -2 at 
4.2 K. The dimensions of the dots defined by the gate pattern are 320 x 320 nm (dot 1, left) and 280 x 280 nm 2 (dot 2, right) 
(b) Schematic diagram of a double dot defined by a combination of dry etching and metallic gates. The carrier concentration 
and mobility of the ungated 2DEG (100 nm below the surface) at 1.6 K and in the dark are 3 x 10 11 cm~ 2 and 8 x 10 5 
cm 2 /(Vs), respectively. The lithographic distance between the ethched trenches (black lines), w, is typically 0.5 fim. The 
effective width of the channel, w c h, can be tuned by voltages on the in-plane gates (IPGs). The gate electrodes are ~40 nm 
wide (l g ) and are separated by Ilc = 160 nm and Icr = 220 nm. A double quantum dot (dot 1, Dl; dot 2, D2) can be formed 
by applying negative gate voltages to gates GL, GC and GR. A microwave field can be applied to the center gate GC. 






(d) 




- 








-*■ — 



Vf 



GL 

FIG. 7 Experimental gray-scale plots of stability diagrams in a device similar to that shown in Fig. |(b) for increasingly 
negative gate voltage on the middle gate electrode, GC. Dark (light) gray-scale corresponds to large (small) current through 
the double dot. GL is swept between -500 mV and -530 mV, GR between -840 mV and -900 mV, GC = -660 mV (a), -670 mV 
(b), -690 mV (c), -700 mV (d), -710 mV (e) and -720 mV (f). 



The stability diagram is of great value in setting up and characterizing a double quantum dot. Figure [7] illustrates 
the process of the creation of a double dot in a device similar to the one in Fig. ||(b). Starting point is the creation 
of a single large dot formed by the outer tunnel barriers, GL and GR. The measured stability diagram (Fig. |7|(a)) 
resembles Fig. @(c). The successive stability diagrams are measured for increasingly negative voltages on the middle 
gate electrode, GC, thus reducing the coupling between the dots. It is clearly seen that going from Fig. ^|(a) to Fig. 
CT, the stability diagram gradually evolves into the characteristic honeycomb structure. The edges of the honeycomb 
cells are visible due to off-resonance current. At the edges of a honeycomb cell, the electrochemical potential of one 
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of the dots aligns with its neighboring lead (see Fig. |). By a process called co-tunneling ( [Averin and Nazarov, 1992[ ) , 
transport can still take place via an intermediate virtual state. Co-tunneling processes are suppressed by increasing 
the tunnel barriers (i.e. making the gate voltages more negative), as can be seen from Fig. [?]. 
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FIG. 8 (a) Experimental gray-scale plot of a stability diagram in the device of Fig. ^(a) at small bias voltage, V = 15 ^eV. 
Dark (light) gray-scale corresponds to large (small) current through the double dot. The dashed lines indicate the honeycomb 
cells, (b) Region within the black rectangle of (a) at large bias voltage, V = 120 ue V. The triple points have grown into 



triangles and show clear resonant tunneling lines (black stripes), as discussed in section 
accentuated by dashed lines. 



II.B.2 



The shape of the triangles is 



Figure |8|(a) shows a detail of a stability diagram obtained in the device shown in Fig. |6|(a). The edges of the 
honeycomb cells are indicated by dashed lines. The triple points within the black square are well separated, whereas 
the other ones are still grown together. To separate also those points outside, the gate voltage on the middle barrier 
has to be tuned towards more negative values. 



B. Non-linear transport regime 

1. Classical theory 

We assume that the bias voltage is applied to the left lead (/i£ = — \e\V) and that the right lead is grounded (/ir 
— 0). The bias voltage is coupled to the double dot through the capacitance of the left lead, Cl, and hence also 
affects the electrostatic energy of the system. The bias dependence can be accounted for by replacing C gl ( 2 )K ;1 (2) 
with C gl (2)V^ 1 ( 2 ) + C L in)V in Eq. (|j), where C il ( 2 ) is the capacitance of the left lead to dot 1 (see Appendix A). 

The conductance regions at finite bias change from triple-points to triangularly shaped regions (Fig. ^J). The 
conditions — \e\V = fiL > Hx, A*i ^ and ybi > [Ir = determine the boundaries of the triangular regions. The 
dimensions of the triangles SV g \ and SV g 2 (see Fig. ||) are related to the applied bias voltage as follows 

a 1 SV gl = ^\e\5V gl = \eV\ 




where ot\ and «2 are the conversion factors between gate voltage and energy. Combining Eqs. (||), (|h]) and ([hi]), we 
can calculate the values of the total capacitances and mutual capacitance C m . 
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FIG. 9 Region within the dotted square of Fig. |2j(a), corresponding to the 'unit cell' of the double-dot stability diagram, at 
finite bias voltage. The solid lines separate the charge domains. Classically, the regions of the stability diagram where current 
flows, are given by the gray triangles. In the case of one discrete level per dot, as in the schematic pictures, resonant tunneling 
is only possible along the side of the triangle that coincides with the dashed line connecting the original triple-points (• and 
O). However, also in this case inelastic tunneling and co-tunneling still contribute to a finite current within the gray triangles. 



2. Quantized states 



For sufficiently large bias voltages, multiple discrete energy levels can enter the bias window. In this case, not only 
ground states, but also excited states contribute to the conductance. For the illustrative case of two levels per dot, 
the four possible alignments of the electrochemical potentials are shown in Fig. [n]. Note that the electrochemical 
potentials are drawn for the situation where one electron is on the double dot and a second one is tunneling on to 
it. Due to Coulomb blockade, not less than one and not more than a total of two electrons is allo wed on the double 
dot. The labeling of the electrochemical potentials, using the notation introduced in section II. A. 2 is straightforward, 
except for tunneling through the excited state of dot 1 in Figs. 10(c),d. Although the second electron is tunneling into 
the lowest level available in dot 1, this level is only accessible because dot 1 is in an excited state. For that reason, 
we choose the label ^1^(2,0) (instead of /ii.o(2,0)). The successive alignment of ground and excited states leads to 
resonances within the conductance triangles, as sh own in Fig. pd[ The off-res onance conductance in the grey triangles 



is due to inelastic processes (Fujisawa et ai, 1998 
199|). Note that V is so large that the two triang 



van dcr Vaart et ai, 199"5j ) and co-tunneling (Averin and Nazarov 
es partly overlap 



Figure |^(b) shows the triple points within the black rectangle of Fig. ||(a) at finite bias. The triangular regions are 
clearly visible as well as the resonances within the triangles. The growth of the triangular regions with increasing 
bias voltage is illustrated in Fig. [l^. Whereas only the ground state resonance is observed in Fig. [l2](a), multiple 
resonances appear within the triang les of Fig. @b). 



III. RESONANT TUNNELING 



In t his section we discuss resona nt tunneling experiments through the double dot of Fig. ^(a) with discrete energy 
levels (van der Vaart et al., 1995). We show that, under appropriate conditions, the resonance widths are only 
determined by the lifetime of the discrete energy states, independent of the electron temperature in the leads. A small 
asymmetric deviation from the Lorentzian resonance shape is attributed to inelastic tunnel processes. 
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FIG. 10 Schematic diagrams showing the possible alignments of the electrochemical potentials in the case of two levels per dot. 
(a) The first electrochemical potentials to align correspond to the ground states of both dots, Gi and G2. (b) When moving 
down the levels in the right dot, the next states to align are the ground state of the left dot, Gi, and the first excited state of 
the right dot, X2- (c) Shifting the levels of the right dot further down, results in transport through the first exited state of the 
left dot, Xi and the right dot ground state, G2. (d) Finally, the excited states, Xi and X2 align. 



The current-voltage (I-V) curves of the single quantum dots in Fig. |l3| provide two clear signatures for the presence 
of both Coulomb blockade effects and discrete levels. A t low bias voltages, the current through the dot is suppressed 
by the Coulomb blockade ( Kouwenhoven et ai, 1997 ). Increasing the bias voltage lifts the blockade. The current 
shows a stepwise increase: each time when a n additional level en t ers the bias window — \e\V, an extra transport 



Johnson et al, 1992| ; | u en 



cha nnel is opened and the current increases ( Foxman et ai, 1993| ; |Gueret et ai, 19*92 ; 

al., 1992| ). Hence, the voltage spacing of the current steps directly reflects the energy spacing of the levels. For the 
average level spacing, 8, we obtain 81 = 125 fieV for dot 1 (upper inset) and 82 = 225 /xeV for dot 2 (lower inset). The 
difference in these two energies reflects the different lithographic sizes of the two dots (see Fig. ^(a)). Accounting for 
the depletion areas, we estimate that dot 1 has an effective diameter of 240 nm and contains about N% = 90 electrons, 
while dot 2 has an effective diameter of 200 nm and contains roughly N2 = 60 electrons. Using the Fermi energy, 
Ep, at bulk density, we estimate Si « 2Ep/Ni = 150 /ieV and 82 ~ 2Ep/N2 = 230 fieV. This is in good agreement 
wit h the estimates obtained from the I-V curves. From the dimensions of the Coulomb diamonds ( Kouwenhoven ei 
al. 



L997; Pothicr et al, 1992) we obtain the charging energies Ec for adding an electron to the dot : Eqi 
meV (dot 2). 



1.1 meV 



(dot ljand E C2 = 1- 

We focus on the role of the discrete levels and consider the charging energies as constant offsets in the transport 
conditions. Figure |l3| shows an I-V curve of the double dot with all three tunnel barriers set in the weak-tunneling 
regime. The Coulomb blockade suppresses the current through the double dot at low bias voltages. At larger bias 
the current shows sharp resonances. The spacing of the resonances is about 250 /ieV. This is of the same order as the 
level spacing in the single dots. 

The same resonances are seen when we sweep the gate voltage. Figure [l4] shows the current through the double 
dot versus the gate voltage on gate 1, V g i, with V = 280 ^eV. This corresponds to a vertical cut through a stability 
diagram as shown in Fig. |l2|. The current shows three groups of sharp resonances separated by regions of zero current 
with a period AV g i = 9 mV in gate voltage V g i- With only dot 1 formed, we observe Coulomb oscillations as a function 
of V g i with the same period AV^i; each period thus corresponds to a change of one electron in dot 1, while keeping 
the number of electrons on dot 2, N2, constant. AV^i corresponds to the horizontal dimension of the honeycomb unit 
cell, as indicated in Fig. || in section II.A.l. 

When elastic tunnel processes are the dominant transport mechanism, the current through the double dot is 
resonantly enhanced only when two levels in dot 1 and 2 align, as explained in section [LB. 2. Tuning the level alignment 
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FIG. 11 Schematic stability diagram corresponding to the finite-bias diagrams of Fig. nfl. The black solid lines within the gray 
triangles correspond, from bottom to top, to the level alignments shown in Fig. [To|(a)-^d), respectively. 
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FIG. 12 Experimental gray-scale plots of stability diagrams in the non-linear regime, obtained in the device of Fig. ^|(a). Dark 
(light) gray-scale corresponds to large (small) current through the double dot. The bias voltage between source and drain 
contacts is 200 //V (a) and 1 mV (b). 



with V or V g i gives rise to the sharp resonances in Figs. [L3| and [14]. Resonant tunneling through the double dot is 
illustrated in the schematic potential landscape of the double dot in Fig. [ll| This figure shows a few of the levels in dot 
1 (levels 1 to 5) and dot 2 (levels a and /?). The electrostatic potentials ipx and <f 2 are tuned in such a way that transport 
through the double dot is possible only via the charge states {N\,N 2 ) -> {N x + 1,N 2 ) -> {N U N 2 + 1) -> (N U N 2 ). 
The finite bias voltage gives an electron from the left reservoir three choices to tunnel into dot 1: it can tunnel to one 
of the unoccupied levels 3, 4 or 5. This changes the electrostatic potential <f\ by the charging energy Eci (the levels 
are drawn at the positions applicable after an electron has occupied one of them) . When dot 1 relaxes to the ground 
state (the incoming electron occupying level 3), the electron can tunnel via level a to the right reservoir. 

Note that if dot 1 does not immediately relax to the ground state, but remains in an excited state, with an electron 
occupying either level 4 or 5, electron transport through the double dot is temporarily blocked. The electron is 
'trapped' within dot 1. Only after relaxation to the ground state a next tunnel event can occur. If the relaxation 
rate is small on the scale of the tunneling rates through the barriers, the inclusion of levels 4 and 5 within the bias 
window could therefore lead to a decrease of the current through the double dot. On the other hand, in case of fast 
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VOLTAGE (mV) 

FIG. 13 I-V curve of the double dot, showing sharp resonances in the current when two discrete levels align. Upper inset: 
I-V curve of dot 1. Lower inset: I-V curve of dot 2. Both insets show a suppression of the current at low volta ges due to the 
Cou lomb b lockade and a stepwise increase of the current due to the discrete energy spectrum of the dot (from ( van der Vaart 
et 4, 19951))- 
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FIG. 14 Current through the double dot versus gate voltage V g i using a bias voltage V = 280 /iV. Inset: The current through 
t he double dot as a functio n of V g 3 with V = 1 mV, showing that the number of resonances increases with bias voltage (from 
(van der Vaart et al, 1995)). 



relaxation, the enhanced tunnel probability when also the levels 4 and 5 lie within the bias window, could lead to an 
increase in the current. Note that next to intra-dot relaxation, also inelastic tunneling fro m level 4 or 5 to eithe r level 
a or (3 can occur. This process is accompanied by emission of a boson (usually phonons ( Fujisawa et al., 1998| ) ) and 
contributes to the off-resonance current in Fig. \L4. 

The resonances in a particular group in Fig. |14| can be identified with the energy diagram of Fig. [l^. The first 
resonance occurs when level 3 aligns with level a (peak 3-a). This corresponds to the rightmost peak in Fig. |lj. 
Increasing — \e\ifi by making V g \ more negative, brings transport off-resonance until level 2 aligns with a (peak 2-a) 
followed by the third peak 3-/3. Continuing to sweep V g \ increases the energy of level 3 above the electrochemical 
potential /i£ of the left reservoir. This blocks transport and removes an electron from dot 1 permanently The next 
group of resonances is observed when V g \ is changed by one Coulomb oscillation period &V g \ (see Fig. |lj). Note that 
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FIG. 15 Schematic potential landscape of the double quantum dot, where [il and [ir denote the electrochemical potentials of 
the left and right reservoirs and V the bias voltage across the do uble dot. The OD states in dot 1 are denoted by levels 1 to 5 
and in dot 2 by levels a and /3 (from (van der Vaart et al, 1995)). 



the number of resonances decreases in the next two groups. Sweeping V g \ also shifts the levels in dot 2, due to a small 
cross-capacitance between gate 1 and 2. Transport is possible until lev el a is shifted above jjj, (V a i < -470 mV) 



The level spacing is obtained by converting gate voltage to energy (Grabert and Dcvorct, 1992) . This yields an 
energy separation of resonances 2-a and 3-a by 70 ^teV, which is the energy separation of levels 2 and 3. In the same 
way we find for levels a and (3 a separation of 200 /xeV. Both values are in good agreement with the typical values we 
found above. On increasing V, we observe that the number of resonances in a particular group increases. The inset 
to Fig. |l4] shows approximately 11 resonances as V g z is swept. These observations are in agreement with the resonant 
tunneling picture of Fig. I5f when V is larger, more levels can align. 
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FIG. 16 Current through the double dot of Fig. ^(a) versus the voltage on gate 3, V g z for bias voltages, V, between 0.2 mV 
(lower trace) and 1.6 mV (upper trace). The traces have been given an offset proportional to their bias voltage for clarity. 
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Generally, the relaxation rate to the ground state is not necessarily higher than the tunnel rate through the dot. In 
Fig. [l6] the amplitude of the ground state resonance (see lower curve) clearly decreases with increasing bias voltage. 
At the same time new resonances appear, having a larger amplitude than the ground state resonance. This implies 
that in general transport through excited states can play a significant role. 




FIG. 17 Enlarged resonance measured in a second device of identical design using a bias voltage of 400 fieV. The data points 
(black dots) are fit to a Lorentzian line shape (soli d line). For comparison w e plot a thermally broadened resonance with a 
fitted temperature T = 34 mK (dashed line) (from (van der Vaart et at, 1995)). 



When a discrete level in dot 1 is at a distance much larger than the thermal energy, ksT, from the electrochemical 
potential of the left lead, dot 1 acts as low-temperature-pass filter for dot 2, such that only cold electrons contribute 
to the current (Kouwenhoven, 1995). If the energy levels in the dots are separated by more than ksT, the occupation 
of excited states becomes suppressed. This effectively leaves the dot at zero temperature. Hence, the double dot 
geometry allows for an accurate measurement of the intrinsic line width of the discrete levels, which is not averaged 
by the Fermi-Dirac distribution of the electrons in the reservoirs. In other words, an energy resolution better than 
ks T can be obtained. T he line shape of a resonance is lorentzian when only elastic tunneling is important (Nazarov 
et a \., 1993 ; 3toof, 1997). In our geometry (assuming a bias voltage sufficiently large such that electrons must tunnel 
from the left lead to the left dot to enter the system, and must tunnel from the right dot to the right lead to leave it 
again) the current is given by 



= e ^ 3M2 , „ N (15) 



'(AW + 5 + |t ia p(2 + £) 



where AE is the energy difference between two discrete energy levels in the two dots, Ti is the tunnel rate from the 
left lead to dot 1, \ti 2 \ is the modulus of the tunnel coupling between the two dots and T 3 is the tunnel rate from dot 
2 to the right lead. Note that for elastic tunneling the resonance width is only determined by the lifetime of the 0D 
states and independent of temperature. 

Figure [l7| shows a single resonance (black dots). The right-hand side of the peak fits very well with the lorentzian 
line shape of Eq. (|l5|) (solid line), while the left-hand side shows a deviation from the lorentzian fit. The only free 
fit parameter is the full width at half maximum, FWHM = 5 fieV. From the maximum current and the width of 
the resonance we find with Eq. ( |l5| ) a tunnel coupling |i 12 | ~ 0.2 /ieV and a tunneling rate between the right dot 
and the right lead T3 « lOfieV. For compa rison, we have fitt ed the resonance with a thermally broadened resonance 
I(AE) ~ cosh(A£72fcT)- 2 (dashed line) ( |Beenakker, 199l[ ). The top is fit very well for T = 34 mK, but there is 
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a large deviation in both tails of the resonance. On the right-hand side, the deviation can be accounted for by the 
lorentzian broadening. At the left-hand side, the deviation consists of two components. The first one is the same 
lorentzian broadening as observed on the right-hand side. The second one, however, is an asymmetric contribution 
only occurring on this side of the resonance. The asymmetric contribution to the current appears at the side where 
an electron tunnels from a higher to a lower electrochemical potential. Upon reversing the sign of V, we find that 
the asymmetry appears at the other side of the resonance. As shown in Ref. ( Fujisawa et al., 1998| ), this is due to 
inelastic tunnel processes. In such a process, an electron can tunnel inelastically and spontaneously emit its energy 
as a photon or a phonon. 



IV. MAGNETIC FIELD SPECTROSCOPY 



In this section we measure the energy evolution versus magnetic field, B, of energy states near the Fermi energy, 
Ep, in the double quantum dot shown in Fig. [|(a). As a function of B and the voltage on gate 3, V g 3, we observe 
crossings and anti-crossings between Coulomb peaks. The resolution is high enough that avoided crossings in the 
spectrum of a quantum dot can be resolved ( postcrkamp et al, 1998a ). To our knowledge, these are the only existing 



data revealing intra-dot level repulsion in a quantum dot system. 

The experiments are performed in the weak coupling limit, such that mixing between quantum states in one dot 
with states in the other dot or in the leads is negligible (see section |v|). We sweep the gate voltages over small ranges 
and focus on a particular charging transition; i.e. transitions between (A^ + ljA^) and (Ni,N2+l) only. Since we 
discuss only one transition at a time, we can, for simplicity, leave out the Coulomb energies from the discussion and 
concentrate on the alignment of discrete e nergy l evels. 

Using the notation introduced in section [I.A.2| , we label the accompanying electrochemical potentials fii, n (Ni, N 2 ) 
for dot 1 and (J-2,n(Ni, N2) for dot 2 (or simply /ii(2), n )- The condition for tunneling between the lowest possible 
states, i.e. from ground state to ground state, is f^ifl(Ni + 1, N 2 ) = ^2.0(^1, N2 + 1). To tunnel from the first excited 
state of dot 1 to the ground state of dot 2, the condition becomes fii i(Ni + 1, N2) = ^2,0(^1, N2 + 1). The changes 
in the electrochemical potential /ii,n+i — fii, n and fi 2 ,n+i — M2,n are typically ~ 150-200 /j,eV. 

Figure |l^(a) shows a typical set of current traces for different magnetic fields while sweeping V g z- The bias voltage 
V = 1.2 mV is such that several discrete levels in each dot are between the electrochemical potentials of the two leads. 
This is similar as in the inset to Fig. |l4|, but now measured for different magnetic fields. The change in Coulomb peak 
position versus B, AVg 3 eak (AB) , is proportional to the difference in the B-evolution of the electrochemical potentials 
Hi,n(B) and H2,n{B) 



Afi(AB) = [fi^ n (B + AB) - Mi,n(5)] - Wn{B + AB) - n 2 ,n( B )] 
ft 



= aAV„% eak (AB) (16) 



where a is the conversion factor between V g 3 and the electrostatic potential of dot 2. Note that if the states fj,i, n (B) 
and fJ. 2 .n(B) have the same B-dependence, the Coulomb peak position does not change. The energy resolution of 
A/i(AB) is ~ 5 /ieV, corresponding to ksT ~ 50 mK. 

The data in Fig. [l|(a) contain several interesting features. First, we observe crossings between different peaks as 
well as anti-crossings (two are indicated by arrows). Second, pairs of peaks exhibit the same B-dependence. These 
are general features that we observe at several charge transitions (that is for several choices of N\,N 2 )- Independent 
measurements on one of the individual dots also show states evolving in pairs below B ~ 0.5 T. The observed pairing 
and (anti-)crossing of the Coulomb peaks in the double dot experiments can then be explained as shown schematically 
in Fig. [l^(b). Suppose two levels in one dot have an anti-crossing in their B-dependence. Then two paired levels in 
the other dot, having the same independence, both probe this anti-crossing. At the points where two Coulomb peaks 
actually cross, two levels in dot 1 align with two levels in dot 2 simultaneously (though only one electron can tunnel 
at a time, due to Coulomb blockade). 

For the interpretation of the data as schematically given in Fig. |l8|(b), tunneling through the excited state of dot 
1 is a key ingredient. If dot 1 would relax to its ground state much faster than the tunnel rate through the barriers 
after an electron has tunneled onto it via we would only observe the two lower traces in the rightmost diagram 
of Fig. [L8|(b). However, our data suggest a relatively slow relaxation rate between the excited and ground state of 
dot 1 . Recent experiments on (single) quantum dots have shown that indeed relaxation times can be of the order of 



fis or longer when relaxation to the ground state involves electron spin flips (Fujisawa et al., 2001). The condition 
for the relaxation rate from the excited state to the ground state in dot 2 is more subtle. The electron can tunnel 
onto dot 2 via ^2,1 and leave dot 2 either directly or after a relaxation process from /124 to 112,0 has occurred (this 
second possibility requires 112,0 > Mil)- For transport through the double dot, the relaxation rate in dot 2 does not 
necessarily need to be slow as well. However, if the relaxation rate would be too high, the anti-crossing as shown in 
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FIG. 18 (a) Current through the double dot sweeping the voltage on gate 3, V g z, at different magnetic fields, B. The curves 
have been given an offset for clarity. From the leftmost to the rightmost curve, B increases from 300 mT to 600 mT in 3 mT 
increments. The conversion factor a between V g 3 and the electrostatic potential of dot 2, a — -63 /xeV/mV, is determined 
through independent measurements from which we deduced the energy scale, indicated by the vertical arrow in the lower left 
corner (a does not change in this magnetic field range), (b) The first two diagrams show how levels may evolve in each of the 
two dots as a function of B. When these four levels are scanned along each other by sweeping V g 3, this results in peak positions 
as sketched in the rightmost diagram. 



the middle diagram of Fig. [L8|(b) would smear out. The constant level spacing — /ii.o in dot 1 (see left diagram of 
Fig. fl8|(b)) could be explained by an exchange energy, e.g. when the upper level would correspond to a spin singlet 
state and the lower level to a spin triplet (Farucha et al., 2000). 



V. MICROWAVE SPECTROSCOPY 



I n this section we present microwave (0-5 GHz) spectroscopy experiments (Fujisawa and Tarucha, 1997a ,g; Kouwcn- 
hovqn et al., 1997; Oostcrkamp et al, 1997 ) on double quantum dots for different coupling and microwave power regimes 



( poster kamp et at, 1998b ; van der Wiel et al., 1999). We use photon assisted tunneling (PAT) processes, as described 
in sections |V.B| and V.C, to measure the energy differences between states in the two dots of the devices shown in 
(b). Depending on the str ength of the inter-dot c oupling, the two dots can f orm ionic-like (Blick et 

Livermore et al., 1996|; |van der Vaart et at, 1995) or covalent-like bonds 



Fig s | 6|(a ) and 

al., 1996 ; Fujisawa and Tarucha, 1997a , b[ 



(Blick et ah, 1998; Oostcrkamp et al., 19981 f fvan der Wiel et al, 1999 ). In the first case, the electrons are localized 
on individual dots, while in the second case, the electrons are delocalized over both dots. The covalent binding leads 
to a symmetric and anti-symmetric state, whose energy difference is proportional to the tunneling strength between 
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the dots. 

For the microwave experiments we make use of a coaxial cable. From room temperature to the lK-pot, a 0.085 
inch semi-rigid Be-Cu (inner and outer conductor) coaxial cable is used. From the lK-pot to the mixing chamber, we 
use a 0.085 inch semi-rigid stainless steel (inner and outer conductor) coax. From the mixing chamber to the sample, 
various types of low attenuation semi-rigid or flexible coaxial cable can be used, since here the thermal conductivity 
is no longer a constraint. Finally, the coaxial cable is capacitively coupled (typically through a 10 pF capacitor) to 
one of the gate electrodes of the sample, usually the center gate. This gate is capacitively coupled to both dots, and 
hence part of the incident power can generate a microwave oscillating potential across the center barrier. 

A. Two-level systems 

So far, we assumed a purely electrostatic coupling between both dots, whereas tunnel coupling was neglected. 
However, when electrons can tunnel coherently from one dot to the other at appreciable rates, the eigenstates become 
delocalized, extending over the entire double dot system. In principle, these are quantum-mechanical many-body 
states of the two coupled dots. It is very difficult to give a full description of such a many-body system. Therefore, 
we discuss here the elementary case of a quantum-mechanical two-level system, which is quite adequate in grasping 
the physics of a tunnel- coupled double dot. Basically, we only take into account the topmost occupied level in each 
dot and neglect the interaction with electrons in lower energy levels. 



We consider a doub le dot consisting of two well-separated dots, described by a total Hamiltonian Ho (Cohen 



Tanioudji et ah, 1977), with eigenstates |0i) and 102) , and eigenenergies E\ and E 2 (Fig. |l9|(a)) 



H |0 X > = E x \4>x) 

Ho |0 2 ) = E 2 |0 2 > (17) 



We introduce a finite tunnel coupling between the levels in both dots described by the Hermitian matrix T, which for 



simplicity ( Cohen- Tannoudji et ai, 1977) we assume to be purely non-diagonal 

f ° *J 2 ) , h 2 =t* 21 , t 21 = \t 21 \e^ (18) 

One obtains a new Hamiltonian, H = Ho + T, with delocalized eigenstates \ips) (symmetric state) and \ip A ) (anti- 
symmetric state), and eigenvalues E$ and E A 

U\^ S ) = E S \^ S ) 

H \x/, A ) = E A • (19) 

The new eigenvalues can be expressed in terms of the eigenvalues of the uncoupled double dot and the tunnel matrix 
elements as follows 



£'s Em - A /J(A£)2 + |t 12 |2 



E A = E M + \j\{^Ef + M 2 (20) 

where Em = + E 2 ) and AE = E\ — E 2 and \t\ 2 \ = \t 2 i\. The eigenstates an d \i>A) m the basis of |0i) and 
|02) are written 

\il> s ) = - sin V^/ 2 100 +cos^/ 2 |0 2 ) 

|V»A> = cosV^/ 2 |0 1 )+sin^/ 2 |0 2 ) (21) 
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(a) 




FIG. 19 Schematic diagrams of a two-level system, (a) Unperturbed energy levels E\ and E2 (solid lines), and energy levels 
belonging to the symmetric state, Es, and the anti-symmetric state, Ea- (b) Energies Es and Ea versus the energy difference 
AE = E\ — E%. For vanishing coupling (|ii2| ~ 0), the levels cross at the origin (dashed straight lines). For non-zero coupling, 
an 'anti-crossing' occurs: the curves belonging to Es and Ea as function of AE are branches of a hyperbola (solid lines) whose 
asymptotes are the unperturbed levels (see also ( Cohen- Tannoudji et ai, 1977)). 



with tan# = 2\ti2\/AE. Figure |l9|(b) shows the eigenenergies of the coupled two-level system as function of AE. The 
renormalized energy difference, AE*, is given by 



AE* =E A -E S = V(A£) 2 + (2M) 2 



(22) 



Note that the effect of the coupling is stronger for small AE, i.e. close to the crossing of the unperturbed energies 
Ei and E2- Where E\ and E2 cross (Ai? = 0), we have an anti-crossing of Ea and Es, with Ea — Es = 2|£i2|, 
the minimum bonding-anti-bonding energy difference. For large AE, the eigenenergies of the coupled double dot 
approach the eigenenergies of the uncoupled dots, E\ and E 2 - The general solution of the time- dependent Schrodinger 
equation can be written in the form 



(23) 



With Eq. (pi]), \if>(t)) can be expressed in terms of \<pi) and l^)- Since \<pi) and 1^2) are not eigenstates of the total 
Hamiltonian H, they are no longer stationary states. If the system is in state \<pi) at time t — (jV'(O)) = \<fii)) the 
probability P\2{t) of finding it in the state \<p2) a t time t is 



PMt) = \(h\m)f 



4|ii 



4|ii 2 | 2 



(AE)- 



sin 



2h 



(24) 
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Equation (EM) describes a coherent charge oscillation in the double dot system. 



B. Photon assisted tunneling in weakly coupled dots 



If the inter-dot coupling is weak, electrons are strongly localized on the individual dots. In section II. A. 2 we saw 
that we expect a resonant current through the double dot system if ul > ui — (12 > I^r- If we represent a weakly 
coupled double dot by a two-level system, we need the discrete energy levels Ei and E2 to align within the bias 
window. We will only consider the discrete, quantum contribution to the electrochemical potentials and therefore 
simply use the discrete level notation E\ and E2 instead of /Ltj and fi2- 

An additional time-varying potential V r oc cos(27r/t) can induce inelastic tunnel events when electrons exchange 
photons of energy hf with the oscillating field (frequency, /, is typically 1-75 GHz in our experiments) . This inelastic 
tun neling wit h discrete energy exchange is known as photon assisted tunneling (PAT) (Blick et ai, 1995; Kouwcnhoven 



et a 



1994a|Jb|). PAT through a single quantum dot with well resolved discrete OD-states is reviewed in Ref. ( |van| 
der |Wiel et ai, 2002| ). PAT is an invaluable spectroscopi c tool for studying the energy spectra of quant um dots. A 



theoretical study of PAT in double dots is given in Refs. ( Hazelzet et ai, 2001 ; Btoof and Nazarov, 1996). A volta ge 



drop VJj C cos(27r/i) across a tunnel barrier modifies the tunnel rate through the barrier as (Tien and Gordon, 1963) 



T(E)= J 2 n(a)T(E + nhf) 



(25) 



Here n = 0, ±1,±2, and T(E) and T(E) are the tunnel rates at energy E with and without an ac voltage, 
respectively. J^(a) is the square of the nth order Bessel function of the first kind, evaluated at a = eV ac /hf, which 
describes the probability that an electron absorbs (n > 0) or emits (n < 0) n photons of energy hf. Thus, the effect 
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FIG. 20 Squared Bessel functions of the first kind, J^(a), for n = 0, ± 1, ± 2. The inset schematically shows the development of 
sidebands of the original energy as a consequence of the microwave field. A positive (negative) n corresponds to the absorption 
(emission) of n photons during the tunnel process. Elastic tunneling corresponds to n — 0. 



of the interaction between a single-electron state with a classical, oscillating field is that the energy state E is split in 
a set of states E + nhf (see inset to Fig. |o|). The power of PAT as a spectroscopic tool lies in the fact that PAT can 
only take place if the energy difference AE equals an integer number times the photon energy hf: AE = nhf, see 
Fig. |2l]. For the multiple photon processes (|n| > 1) to take place, the microwave power needs to be sufficiently large. 

To use PAT as a spectroscopic tool, we can make use of the configurations shown in Fig. ^l|. In the pumping 
configuration (Brunc et ai, 1997; Stafford and Wingreen, 1996| ), the double dot is operated at zero bias voltage. 
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Absorption of a photon with energy hf = AE leads to pumping of an electron from left to right (Fig. |2l](a)) or 
vice versa (Fig. |2lj(b)). The advantage of this configuration is that relaxation due to spontaneous emission does not 
contribute to the current. Figure |22| schematically shows how the honeycomb unit cell of Fig. || changes in the presence 
of a microwave field. 

Alternatively, the double dot can be operated in the large bias regime as depicted in F ig. pl|(c),d. In this regim e, 
in the case of weak coupling with \ti2 \ <C AE, hf, KTl.r, the dc PAT current is given by ( Stoof and Nazarov, 199(i| ) 
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FIG. 21 Schematic electrochemical potential diagrams of photon assisted tunneling (PAT) in a weakly coupled double quantum 
dot. The upper diagrams (a) and (b) show absorption of a photon with energy hf = AE in the so called pumping configuration. 
Although V — 0, an electron can tunnel from left to right through the dot (a), or vice versa (b). The lower diagrams show 
absorption (c) and (stimulated) emission (d) of a photon with energy hf — AE in the large bias configuration. 



IpAT = e\t 12 \ 2 T R 3 2 n (a)/(-rl+(n27rf-AE/h) 2 ) (26) 

n— — oo 

The current is composed of a number of satellite peaks, separated by the photon energy hf and all with width T^. 
Note that the satellite peaks can become of the same order of magnitude as the main resonance, but that they have 
a smaller width than the main resonance. The PAT experiment described below, is performed on a weakly coupled 
dot in the large bias regime. 

The double dot, shown in Fig. |^(a), is tuned such that only one level in each dot contributes to electron transport. 
The gate voltages are used to shift the level in dot 1 and in dot 2. The resonance in the lowest trace in Fig. ^3|(a) arises 
from the alignment of the two levels. The other traces are measured while applying a microwave signal. The satellite 
resonances are due to PAT processes which involve the emission (left satellite peak) or absorption (right satellite peak) 
of one photon. Figure |2^(b) shows that the energy separation of the satellite peaks from the main peak, AE, depends 
linearly on frequency between 1 and 50 GHz. As we will discuss below, this linearity implies that the tunnel coupling 
is negligible. The electrons are localized on the individual dots and they have an ionic bonding. The line proportional 
to 2hf is taken from data at higher microwave powers where electrons absorb or emit two photons during tunneling. 

As the microwave power is increased, more satellite peaks appear corresponding to the absorption of multiple 
photons, which are observed up to n = 11 (see Fig. p4| ). A high power microwave field strongly perturbs tunneling. 
This is reflected by the non-linear dependence of the peak heights on microwave power. In the right inset to Fig. 
pi] the peak heights of the main peak and the first four photon satellite peaks are shown, which agree well with the 
expected squared Bessel function behavior shown in Fig. EG. 
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FIG. 22 Schematic of a double quantum dot stability diagram in the weak coupling and linear transport regime irradiated by 
microwaves with frequency /. Next to the triple points, finite current is measured along the PAT lines at a distance Spat 
from the line 0) = /^(O, 1), where A_E = hf. The various PAT processes are illustrated by the electrochemical potential 

diagrams. 



C. Photon assisted tunneling in strongly coupled dots 

The large bias configuration of Figs. |2l](c),d was successfully employed to study PAT in a weakly coupled double 
dot system. For the microwave spectroscopy of a strongly coupled double dot we will make use of the pumping 
configuration shown in Fig. ^5|. With increasing the coupling between the dots, the spontaneous emission rate from 
the higher level to the lower one increases as well. The advantage of the pumping configuration is that these processes 
are 'filtered out' and do not contribute to the current. 

Figure ^5|(a) schematically shows the symmetric and anti-symmetric states in the double dot for AE = 0. When 
microwave radiation is applied with a frequency such that hf — AE* — Ea — Eg = 2|ti2|, electrons are pumped 
from the left lead to the right lead and vice versa. Since the weight of the symmetric and anti-symmetric wave 
function is distributed equally over both dots (see Fig. |25|(a)), there is no net current. However, if we detune the 
levels (|AS| > 0) the weight gets distributed asymmetrically as shown in Figs. ^5|(b),c and a net current is generated 
by applying microwave radiation matching hf = AE*. Note that for frequencies hf < 2 ]"ix2 1 no PAT is possible. 

Multi-photon processes occur if the condition AE* = nhf (|n| > 1) is met. Besides allowing for these higher order 
photon processes, a high power microwave field also renormalizes the tunnel coupling to a smaller value. The energy 
splitting AE* now becomes 



AE* = V(A£) 2 + [2J (a)|ii 2 |] 2 . (27) 
The experiments for strong inter-dot coupling are performed on the device shown in Fig. ^(b). To single out 



the current only due to microwaves, we operate the device as an electron pump driven by photons ( Brunc et al. 
199|; Stafford and~~WTn" green, 1996| ) (see the diagrams in Fig. p5[ ). An electron is excited from the bonding to the 



anti-bonding state if the condition hf = AE* is fulfilled, or conversely 

AE = v/^/) 2 -(2Jo(a)M) 2 . (28) 

Figure ^6| schematically shows the stability diagram for a strongly coupled double dot in the presence of a microwave 
field. Here we assume that the symmetric and anti-symmetric states can only be occupied by a single electron. In 
other words, we assume spinless electrons. The triple points of the weakly coupled double dot, denoted by • and O, 
develop into a black and a white crescent, respectively. The length of these crescents increases with | i x 2 1 - 
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FIG. 23 Weakly coupled double quantum dot in the low microwave power regime, (a) The upper schematic pictures illustrate 
three configurations of the discrete energy level in the left dot relative to the level in the right dot (thick solid lines). The 
electrochemical potentials of the leads are indicated by thinner solid lines. The bottom curve shows the current as a function 
of the voltage on gate 1, V g i, (see Fig. |^(b)) for source-drain voltage, V = 500 fiV without applying microwaves. A single 
resonance occurs when two levels align. The other curves, which have been offset for clarity, show the current when microwaves 
with frequency / from 4 to 10 GHz are applied. Now, two additional satellite resonances occur when the two levels are exactly 
a photon energy apart. The corresponding photon-assisted tunneling processes are illustrated in the upper diagrams, (b) 
Distance between main resonance and first two satellites as a function of the applied frequency from 1 to 50 GHz. The distance 
is transferred to energy through AE = kAV 9 i where k is the appropriate capacitance ratio for our device that converts gate 
voltage to energy. The agreement between data points and the two solid lines, which have slopes of h and 2h, demonstrates 
that we observe the expected linear frequency dependence of the one and two photon processes. 



Moving along the dotted line connecting the crescents from lower left to upper right, first the symmetric state aligns 
with the electrochemical potentials of the leads (at the black crescent). Current through the double dot is possible 
via the electron transfer process of Fig. ||(d). At the black square, the anti-symmetric state aligns with the leads. 
However, current is blocked, since an extra electron is already added to the double dot and the charging energy Ec m is 
not available yet. At the white square, the electrochemical potential for adding the second electron to the symmetric 
state aligns with the leads. As we assumed single occupation of the delocalized states, current is blocked here as well. 
When arriving at the white crescent, the electrochemical potential for adding the second electron to the double dot 
in the (empty) anti-symmetric state becomes available. This enables the hole transfer process of Fig. ||(d). 

The black bars at a distance Spat from the dotted line connecting the crescents, denote the places where AE* = hf 
and PAT occurs. Note that PAT is only possible if the photon energy exceeds the coupling energy, hf > 2 1 *i2 1 - 
Depending on the level configuration, pumping results in a negative or positive contribution to the current, as shown 
in Fig. ^7](a) (here we choose I > for an electron moving from left to right). Figure H?](b) shows a part of the 
corresponding stability diagram between two triple points, clearly showing the energy regions of constant charge and 
extra transport lines due to PAT. 

Figure ^8] shows measured current traces as a function of the uncoupled energy splitting AE, where from top to 
bottom the applied microwave frequency is decreased from 17 to 7.5 GHz in 0.5 GHz steps. At the highest frequencies, 
the distance between the pumping peaks is close to 2AE. However, the peak distance decreases faster than linearly 
as the frequency is lowered; in fact the peaks follow the dotted hyperbola rather than the dashed straight lines. The 
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FIG. 24 Weakly coupled double quantum dot in the high microwave power regime. The main graph shows current versus gate 
voltage. The dashed curve is without microwaves and contains only the main resonance. The solid curves are taken at 8 GHz 
for increasing microwave powers resulting in an increasing number of satellite peaks. At the right side of the main peak, these 
correspond to photon absorption. V = 700 yN and the photon energy hf = 32 /ieV at 8 GHz. At the highest power we observe 
11 satellite peaks, demonstrating multiple photon absorption, (left inset) Schematic diagram showing multi-photon absorption, 
(right inset) Height of the first four satellite peaks as a function of the microwave amplitude. The observed dependence shows 
the expected Bessel function behavior given in Fig. EB. 




FIG. 25 Schematic electrochemical potential diagrams of PAT in a strongly coupled double quantum dot in the pumping 
configuration. The diagrams show the symmetric state with wave function if>s and energy Es (lower dashed line), and the 
anti-symmetric state with wave function tpA and energy Ea (upper dashed line) in combination with the eigenenergies E\ and 
E2 of the weakly coupled double dot (solid lines), (a) Ei = E2 — 0, AE* = Ea — Es = 2[ti2|. Irradiation with photons 
with energy hf — 2|ti2| leads to PAT, but the net current through the double dot remains zero, (b) By lowering — \e\tpi and 
increasing — \e\(p% the weight of the wave functions is redistributed such that net electron transport from left to right occurs, 
(c) By increasing — \e\<p\ and lowering — \e\1p2 the weight of the wave functions is redistributed such that net electron transport 
from right to left occurs. 



distance goes to zero when the frequency approaches the minimum energy gap between bonding and anti-bonding 
states, hf = 2|ti 2 |. The coupling between the dots can be decreased by changing the gate voltage on the center gate 
to more negative values, or by applying a magnetic field perpendicular to the sample. In Fig. py(a) we plot half the 
spacing between the positive and negative satellite peaks as a function of frequency. The microwave power is kept as 
low as possible in order to meet the condition eV ac <C hf. In that case J 2 , (a) ~ 1 and the general relation Eq. ( pgj ) 
reduces to 



AE = V(ft/) 2 -(2|ti 2 |) 2 . (29) 
Different symbols correspond to different center gate voltage settings and magnetic fields. The solid lines are fits 
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FIG. 26 Schematic of a double quantum dot stability diagram in the strong coupling and linear transport regime irradiated 
by microwaves with frequency /. The symmetric and anti-symmetric states are assumed to be singly occupied. At the solid 
hyperbolic lines either fis or [ia equals zero, marking the separation of charge domains. At the dashed hyperbolic lines either 
fis or ha equals zero as well, but electron transport is blocked. The triple points of the weakly coupled double dot (• and 
O) develop in to the black and white crescent, respectively. At the position of the black and white square no current occurs, 
as explained in the text. At a distance Spat from the dotted line connecting the crescents, AE* — hf and PAT occurs. The 
various configurations of the electrochemical potentials are also illustrated. 



to Eq. (|29|). The good agreement with Eq. ( |29| ) demonstrates the control over the formation of a covalent bonding 
between the two dots and that the condition eV ac <C hf is satisfied. 

We now discuss t he case eV ac > hf. As ca n be seen in Fig. Jg(a) deviates from 1 in this case and cannot be 
neglected as before ( van der Wiel et al., 199S ). In Fig. ||^(b) we show the power dependence for the case of a coupling 
of 60 jtxeV and a microwave frequency of 16 GHz, as indicated by the circle in Fig. ^(a) (similar results have been 
obtained for other couplings and microwave frequencies). The inset to Fig. ^(b) shows the measured PAT current 
as a function of AE for different powers. The absolute value of the microwave power at the position of the double 
quantum dot is unknown. Therefore, we use a relative microwave power scale, which is expressed in terms of the 
attenuation of the microwave source signal. The positions of the PAT peaks at the lowest power are indicated with 
two dashed lines. Increasing the microwave power from the lowest value, the PAT peak separation becomes larger, 
which is in agreement with Eq. (Eq). For higher powers, multi-photon processes can also take place, which result in 
extra current peaks. Figure E9(bjshows half the PAT peak separation energy as function of the relative microwave 
power. The solid line is a fit with Eq. (p8|), / = 16 GHz, 2|tia| = 60 fieV. Because of the relative power scale, the 
fitting curve has been adjusted horizontally to give the best fit. We thus see that the microwave power effectively 
reduces the coupling between the dots. This is further illustrated by the vertical dotted line in Fig. p9|(a) at / = 
16 GHz. At -33 dB the energy separation equals hf, which implies that the higher microwave power has effectively 
reduced the coupling between the dots. We can obtain an estimate of the power by noting that 3q(o) has its first 
zero for a = eV ac /hf = 2.4 and hence V ac = 0.16 mV. 



VI. CONCLUSIONS 

By coupling two quantum dots in series, we obtain a system with fundamentally different behavior and possibilities 
in comparison to a single quantum dot. In this review we have discussed the superiority of a double quantum dot 
system in determining the intrinsic lifetime of quantum states and in probing intra-dot level repulsion. Next to the 
added value as a spectroscopic instrument, the double dot manifests itself as an artificial molecule. By changing the 
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FIG. 27 Strongly coupled double dot in the low-power regime, (a) Current through the double dot as function of the energy 
difference between the level in the left and the right dot. The current trace is taken from the stability diagram in (b) at the 
position indicated by the arrow. The diagrams depict the discrete levels E\ and E2 in the two dots for the case that the 
coupling is weak (solid lines) and the bonding and anti-bonding states in the case of strong coupling (dotted lines). The PAT 
processes leading to a negative (left diagram) and a positive current (right diagram) are indicated, (b) Gray-scale plot of the 
current through the double dot versus the energy level difference, AE, and the total energy, E t . The bias voltage is 6 fiV and 
the applied microwave frequency is 16 GHz such that hf = 66 ^teV. The dashed lines divide the stability diagram in 4 regions 
of stable electron numbers. In between the two triple points clear features of photon-assisted tunneling are seen. The black 
arrow indicates the position of the trace shown in (a). 



inter-dot coupling, we have been able to tune the double dot from an ionic-like bonded to a covalent-like bonded 
molecule. 

Now that the ability to create and manipulate double quantum dots has been shown, the next challenge lies in 
the study and time-control of coherent phenomena in these syst ems. Double quantum dots have been suggested as 
possible candidates for building blocks of a quantum computer ( Loss and DiVinccnzo, 199S| ). We have shown that 
it is indeed possible to coherently couple dots, and that one can induce transitions between the extended states. 
The next crucial step towards quantum logic gates is to show that the coherence of the superposition is preserved 
on time scales much longer than the time needed for manipulating the electron wave functions. The time-resolved 
measurement of coherent charge oscillations in double quantum dots using pulsed gate voltages, will be an essential 
step in determining the dephasing time in these systems. 

In addition, the role of the electron spin and its relaxation and dephasing times (T\ and T2, respectively) need to 
be characterized, as they are essential in the proposed quantum bit schemes based on coupled electron spins in double 



28 




AE (neV) 

FIG. 28 Measured pumped current through the strongly coupled double dot. Gates 1 and 3 are swept simultaneously in such 
a way that we vary the energy difference AE. The different traces are taken at different microwave frequencies, and are offset 
such that the right vertical axis gives the frequency. The main resonance is absent as we have set V = 0. The satellite peaks 
typically have an amplitude of 0.5 pA. For weakly coupled dots the satellite peaks are expected to move linearly with frequency, 
thereby following the straight dashed lines. In contrast, we observe that the satellite peaks follow the fitted dotted hyperbola 
hf = [AE 2 + (2|£i2 1) 2 ] 1//2 using the coupling \ti2\ as a fitting parameter. 



dots ( Loss and DiVinccnzo, 1998| ). The relaxatio n time T± in qua ntum dots has been experimentally shown to exceed 
10 - 100 /is if the relaxation involves a spin-flip ( Fujisawa, 2002 ). The determination of the spin dephasing time T2, 
which is the most relevant timescale for quantum computing purposes, will be an experimental challenge for the near 
future. We conclude that our work on double dots so far, in combination with the results on the spin relaxation times 
in this kind of systems, forms a promising point of departure for further study on the suitability of double quantum 
dots as quantum coherent devices. 
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FIG. 29 (a) Half the spacing in gate voltage between the positive and negative satellite peaks as a function of frequency for 
eVac <C hf. Gate voltage spacing has been transferred to energy difference AE (see also figure caption Fig. ^3|(b)). Different 
curves correspond to different tunnel couplings | £12 1 - Solid lines are theoretical fits to Eq. (E9I). In the limit of weak coupling, 
this reduces to AE — hf, which is indicated by the dashed line. The resulting values for 2 1 £ 12 1 are given in the figure. The 
coupling is varied by applying different voltages to the center gate or by changing the magnetic field (♦ B — 3.3 T; ■ B — 2.2 
T; other curves B = T). The circle marks a coupling of 60 ^teV and frequency of 16 GHz (dotted line) . (b) Strongly coupled 
double dot (2|tia| = 60 AteV) in the high microwave power regime for / = 16 GHz (along dotted line in (a)). The inset shows 
the measured PAT current as a function of AE for different powers. The positions of the PAT peaks at the lowest power are 
indicated with two dashed lines. The PAT peak separation becomes larger for increasing microwave power. For higher powers, 
multi-photon processes can also take place, which result in extra current peaks. In the main part, half the PAT peak separation 
in energy as function of the relative microwave power is shown. The solid line is a fit to Eq. (^), / = 16 GHz, 2 |i 12 1 = 60 /xeV. 
Because of the relative power scale, the fitting curve has been adjusted horizontally to obtain the best fit. 



APPENDIX: Electrostatic energy of quantum dots 



In this appendix we derive the electrostatic energy of a single and double quantum dot system. Before addressing 
these specific systems, wc briefly discuss the method followed 5 . 



1. Electrostatics of a system of N conductors 

Consider a system consisting of N conductors. A capacitance can be defined between each conductor and every 
other conductor as well as a capacitance from each of the N conductors to ground. This results in a total of N(N+ 1)/2 
capacitors. The capacitor between node j and node k has a capacitance Cjk and stores a charge q^. The total charge 
on node j is the sum of the charges on all of the capacitors connected to node j 

N N 

Qi = X>* = X^fVj - (3°) 

fe=0 fc=0 

Here Vj is the electrostatic potential of node j and ground is defined to be at zero potential, Vq = 0. The charges on 
the nodes are linear functions of the potentials of the nodes so this can be expressed more compactly in matrix form 

^-CV" (31) 



For a discussion of the electrostatics of a charging network see 



http://qt.tn.tudclft.n1/~ 



hadlcy/set/electrostatics.html. 
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where C is called the capacitance matrix. A diagonal element of the capacitance matrix, Cjj , is the total capacitance 
of node j 

N 

k=0,k^j 

An off-diagonal element of the capacitance matrix is minus the capacitance between node j and node k, Cjk — Ckj = 
—Cjk- The electrostatic energy of this system of conductors is the sum of the electrostatic energy stored on the 
N(N + l)/2 capacitors and can be conveniently expressed using the capacitance matrix 

U = l^. C V = ^V-^^^-C- 1 ^ (33) 

Voltage sources can be included in the network by treating them as nodes with large capacitances to ground and 
large charges on them such that V — Q/C. In this case, it is numerically difficult to compute the inverse of the 
capacitance matrix since it contains large elements. However, it is not necessary to invert the entire capacitance 
matrix since the voltages on the voltage sources are already known. Only the voltages on the other nodes need to be 
determined. These voltages can be determined by writing the relation between the charges and the voltages as 



(34) 



Here ~C$ c and V* c are the charges and the voltages on the charge nodes, ~C$ v and V* v are the charges and the voltages 
on the voltage sources, and the capacitance matrix has been expressed in terms of four sub-matrices. The voltages 
on the charge nodes are then 

V* c = C- 1 ^ - Ccv (35) 
and the electrostatic energy can be calculated with Eq. (133). 



2. Single quantum dot 




FIG. 30 Network of capacitors and voltage nodes used to calculate the electrostatic energy of a single quantum dot. 



We write the total charge Qi on the dot as the sum of the charges on all the capacitors connected to the dot (see Fig. 
1) 

Qx = C L {V X - V L ) + C g {V x - V g ) + C fl (Vi - V R ) => 

Qi + C L V L + C g V g + C R V R = CiV x (36) 

where C\ is the total capacitance coupled to the dot, C\ = Cl + C g + Cr. The capacitance matrix C cc only has one 
element. Using Eq. ( |33"| ) and substituting Qi — —(N± — A )|e|, we find 

u{Ni) = l-(Ni - N )\e\ + C L V L + C g V g + C R V R \ 2 (3?) 

where Nq is the number of electrons on the dot when all voltage sources are zero, which compensates the positive 
background charge originating from donors in the heterostructure. 
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3. Double quantum dot 




FIG. 31 Network of capacitors and voltage nodes used to calculate the electrostatic energy of a double quantum dot. 



We write the total charge Qx(2) on dot 1(2) as the sum of the charges on all the capacitors connected to dot 1(2) (see 
Fig. H) 

Qi = C L (Vi - V L ) + C gl {Vi - V gX ) + C m {Vi - V 2 ) 

Q 2 = C R {V 2 - V R ) + C g2 (V 2 - V g2 ) + C m {V 2 - V x ) (38) 



We can write this as 





( Cj 


~ C m \ 




)- 


v ~C m 


c 2 J 





Ql + C L V L + C gl V gl 

Q2 + C R V R + C g2 V g2 
where C 2 = Cr + C g2 + C m . The above expression in the form of Eq. ([35]) reads 



(39) 



Vx\ _ 1 / C 2 C m \ (Qx + C L V L + C gl V g x\ , , 

V 2 J CxC 2 - CI {C m Cx J \Q 2 + CrVr + C g2 V g2 ) [W> 

The electrostatic energy of the double dot system can now be calculated using Eq. ([33]). For the case Vl — Vr — 
and Qi(2) = — iVi(2)|e| this becomes 

U(Nx,N 2 ) = ^NfE C x + l -N\E C2 + NxN 2 E Cm + f(V gl , V g2 ) (41) 
f(V g x,V g2 ) = -^-{CgxVgx{NxE C x + N 2 E Cm ) + C g2 V g2 (NxE Cm + N 2 E C2 )} 



- e 



with 



^{2 C 9^ V 9^ Eci + 2 C g2 V g2 E C2 + C g xV g xC g2 V g2 E Cm } 



1 C 2 2 Cx 2 C v 



Eci - e — — — T ; E C2 - e ——^ — r \ E Cm - e — — — T (42) 

OiC/2 — O m O1O2 — O m O1O2 — o m 
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